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G\ ! Abstract 
O 
O 

Considering four-neutrino schemes of type 3+1, we identify four small re- 
gions of the neutrino mixing parameter space compatible with all data. As- 
suming a small mixing between the sterile neutrino and the isolated mass 



43 

a: 



eigenstate we show that large Un — ► v T and v e — ► v T transitions are predicted 
in short-baseline experiments and could be observed in the near future in ded- 



[ icated experiments. We discuss also implications for solar, atmospheric and 

long-baseline neutrino experiments and we present a formalism that allows to 
describe in 3+1 schemes atmospheric neutrino oscillations, long-baseline v„ 
^ . disappearance and — * v T transitions in matter. 

c3 ' PACS numbers: 14.60.St 
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I. INTRODUCTION 



Neutrino oscillation experiments are powerful probes of neutrino masses (see [ITj— . At 
present three types of neutrino oscillation experiments have obtained positive results: solar 
and atmospheric experiments and one short-baseline experiment. 

All solar neutrino experiments (Homestake J7|, Kamiokande ||, GALLEX ||, SAGE 
T(l , Super-Kamiokande [] l]|) have found a deficit in the flux of electron neutrinos on Earth 



with respect to the Standard Solar Model prediction [|12|], which constitutes an indication 
in favor of oscillations of electron neutrinos into other states. Although no direct proof of 
these oscillations exist at present, there is a convincing evidence that the deficit of solar z/ e 's 
is due to neutrino physics, of which neutrino oscillations is the simplest and most natural 
phenomenon. Hopefully, the issue will be definitively settled in a few years by the new 
generation of neutrino oscillation experiments (GNO SNO Jl4]], Borexino ICARUS 
]T§ and others |T7[). 

Several atmospheric neutrino experiments (Kamiokande flTB| , 1MB [HI], Super- 
Kamiokande [|g , Soudan 2 0, MACRO 
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(Kamiokande 0, 1MB 
have found an anomalous ratio of the events 



generated by muon and electron neutrinos and an anomalous angular dependence of the 
events generated by muon neutrinos. Although so far the oscillation pattern has not been 
observed, this is considered as an evidence in favor of oscillations of muon neutrinos into tau 
or sterile neutrinos. Transitions of muon neutrinos into electron neutrinos are strongly disfa- 
vored by the bounds established by the long-baseline v e disappearance experiments CHOOZ 

m 



E3I and Palo Verde 



and by the fact that Super-Kamiokande data do not show any 
anomalous angular dependence of the events generated by atmospheric electron neutrinos. 
The present data of the Super-Kamiokande [20] and MACRO [^2 disfavor pure transitions 
of atmospheric muon neutrinos into sterile states. There are good chances that the issue will 
be clarified in a definite way in the near future by long-baseline experiments with muon neu- 

MINOS fH, OPERA |7|, ICARUS pf) and by new atmospheric 
and others p9[). The K2K experiment have already 



trino beams (K2K || 
neutrino experiments (MONOLITH 
obtained some preliminary indication of a possible transition of muon neutrinos into other 
states EOT . 



The third evidence in favor of neutrino oscillations has been found in the short-baseline 
LSND experiment where an excess of e + events have been observed. If interpreted in 
terms of — » v e neutrino oscillations, this excess corresponds to an oscillation probability 
of (2.5 ± 0.6 ± 0.4) x 10~ 3 . This result has not been confirmed by other experiments (many 
other less sensitive short-baseline experiments have not found any signal of neutrino oscilla- 



tions; the strongest bounds have been obtained in the Bugey p2| , CDHS p3[ , CCFR pi 
BNL-E776 [0, KARMEN GJ, CHORUS |J7|, NOMAD J5§ experiments), but it is very 



interesting and important, because it is the only existing evidence of neutrino flavor tran- 
sition from one state to another and because such transitions could be explored with high 
accuracy in future short-baseline experiments. In a few years the MiniBooNE experiment 
will check the LSND signal in terms of neutrino oscillations. 

H0|,|41I1) that the three indications in favor of neutrino oscillations 



It is well known (see 



need at least three different neutrino mass-squared differences (Am s), arranged in the 



hierarchical order Am| UN <C Am\ TM <c Amg BL 



This means that at least four massive 
neutrinos must exist. Here we consider the minimal case of four massive neutrinos, that 
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has been considered recently by many authors (see [^-|52|,|6|,|40[| and references therein), 



whose flavor basis is constituted by the standard active neutrinos v e 



and a sterile 



neutrino v s . Figure [1] shows the six possible four-neutrino schemes that can accommodate 
the observed hierarchy of Am 2 's. These six schemes are divided in two classes: 3+1 and 2+2. 
In the 3+1 schemes there is a group of three neutrino masses separated from an isolated 
mass by the LSND gap of the order of 1 eV, such that the largest mass-squared difference, 
Am^ (where Amjy = m\ — m 2 and rrik, rrij are neutrino masses, with k,j = 1,...,4), 
generates the oscillations observed in the LSND experiment. In the 2+2 schemes there 
are two couples of close mass eigenstates separated by the LSND gap. The numbering of 
the mass eigenvalues in Fig. |1] is conveniently chosen in order to have always solar neutrino 

and short-baseline (SBL) oscillations generated by 



oscillations generated by Am^ 



^ m SUN 



|Am^| 
jAmfJ 



|Am| 2 | 
|Am 2 2 | 



Am SBL 



[we have also [Am^-J 



|a?t4| 



in 3+1 schemes and |Am 2 



in 



in 2+2 schemes). In 3+1 schemes atmospheric neutrino oscillations are 
generated by lAmU ~ |Amg 2 | = Am\ TM , whereas in 2+2 schemes they are generated by 



|Am 2 3 | 



^ m ATM' 



It has been shown that in the framework of four-neutrino mixing the 3+1 schemes are 



disfavored by the experimental data, with respect to the 2+2 schemes [^2]- |4^j40|1 . However, 
in a recent paper Barger, Kayser, Learned, Weiler and Whisnant [52| noticed that the new 
99% CL allowed region obtained recently in the LSND experiment and presented at the 
Neutrino 2000 conference allows the existence 3+1 schemes (see also f49fl ) in four small 
regions of the amplitude A^ e of — > v e and — > v e oscillations (A^ is equivalent to the 
usual sin 2 2$ in the two-generation case). These regions, to be derived later (Eqs. (|3.1| )- (|3.6| )) 
are shown in Fig. 0, enclosed by very thick solid lines. They lay at 



Rl: 


Am 2 ^ 


~ 0.25 eV 2 


R2: 


lAm 2 ^ 


~ 0.9 eV 2 , 


R3: 


Am 2 ^ 


~ 1.7eV 2 , 


R4: 


\Am 2 41 \ 


~ 6 eV 2 . 



Barger, Kayser, Learned, Weiler and Whisnant explored the phenomenological conse- 
quences of the assumption 



U„A 2 < 1 



;i-2) 



(here U a )~ are the elements of the 4x4 neutrino mixing matrix with a = s,e,fi,r and 
k = 1,...,4). Such a scheme is attractive because it represents a perturbation of the 
standard three-neutrino mixing in which a mass eigenstate is added, that mixes mainly 
with the new sterile neutrino v s and very weakly with the standard active neutrinos u e , 
and v T . In this case, the usual phenomenology of three-neutrino mixing in solar and 
atmospheric neutrino oscillation experiments is practically unchanged. The atmospheric 
neutrino anomaly would be explained by dominant — > v T transitions, with possible sub- 
dominant i/p ^ v e transitions constrained by the CHOOZ bound |53|j54] , |50| . The solar 
neutrino problem would be explained by an approximately equal mixture of v e — > and 
v e — > v T transitions p4|^5 . 
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Here, we consider another possibility that, as we will see, predicts relative large — > v T 
and v e — > v T transitions in short-baseline neutrino oscillation experiments, that could be 
observed in the near future. We consider the 3+1 schemes with 

|C/ s4 | 2 <l. (1.3) 

This could be obtained, for example, in the hierarchical scheme I (see Fig. [I]) with an 
appropriate symmetry keeping the sterile neutrino very light, i.e. mostly mixed with the 
lightest mass eigenstates. Notice that nothing forbids |t/ S 4| 2 to be even zero exactly. 



II. 3+1 SCHEMES WITH \U s4 \ 2 « 1 

From the assumption ( |1.3| ), since the amplitude of v a — > v$ and vp — > v a oscillations 
in short-baseline neutrino oscillation experiments (equivalent to the usual sin 2 2$ in the 
two-generation case) is given by [[|^.[J 

A aj 3 = A Pa = 4\U a4 \ 2 \Up 4 \ 2 (a, (3 = s,e,/x,r) , (2.1) 

we have 

A as <^l (a = e,fi,r), (2.2) 

i. e. the transitions from active to sterile neutrinos in short-baseline experiments are strongly 
suppressed. In the following we will neglect them. 

Let us consider now the oscillation amplitude (equivalent to the usual two-generation 
sin 2 2$) in short-baseline v a disappearance experiments, which is given by [42,|] 



B a = \U a4 \ 2 (1 - \U a4 \ 2 ) . (2.3) 

In general the oscillation amplitude in v a disappearance experiments is related to the am- 
plitude of v a — > vp oscillations by the relation 

B a = J2 A ^, (2-4) 
that quantify the conservation of probability. Using Eq. ( j2.2|) , Eq. (|2.4|) gives 



A eT ~B e - . (2.5) 

A^r -B^- A^ e . (2.6) 

Therefore, we have 

B e ~ A^ £ A eT ^B e - A^ e . (2.7) 

nmin /imax <" a nmax /imin /o o\ 

/i ~ A iie £ a ht r^i By. - Ay £ . (2.8) 

Let us determine Bf n , B™ ax , B™ in , 5™ ax , A™ in , A™ x from the results of short-baseline 
experiments. 
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III. GENERAL BOUNDS IN 3+1 SCHEMES AND 

TRANSITIONS 



SHORT-BASELINE 



The values of B 
experiments (notice that B 

[Am^-J in the LSND-allowed region are given by the exclusion curves of the Bugey JHJ and 
CHOOZ p3) reactor v e disappearance experiments and the exclusion curve of the CDHS 



" HX and _B™ ax are given by the exclusion plots of v e and disappearance 
luax and _B™ ax depend on [Am 2 ^). The most stringent bounds for 



accelerator disappearance experiment |33 



The bounds B e < Bf™ and B^ 



; , together with the results of solar and atmo- 
spheric neutrino experiments imply that \U e ^\ 2 and [t^] 2 are small [42], 43], 40]]: 



f 1 



u e ,\ 2 < \u, 



e4 1 max 



and 



IM 2 < 



\U, 



/"4 1 max ' 



where \U e4 \ max and I I max are g iven by 

1 



u, 



2 

e4 1 max 



\u, 



A*4 1 max 



2 
min 



1 - Bl 



1 Dmax 

ft 



0.55 



(3.1) 

(3.2) 
(3.3) 



The number 0.55 comes |44[ from the up-down asymmetry of multi-GeV muon-like events 



measured in the Super-Kamiokande experiment pD| . The values of |£44| max an d W, 
functions of jAm^l are shown, respectively, in Figs. ^| and |j. 
In the 3+1 schemes we have 

A IU . = 4|[/ /Lt4 | 2 |C/ ( 12 



M 4| 



as 



'■fie 



e4| j 



and the bounds (|3.1| ) imply that [42 



A^ e < A\U, 



A*4l 



el 



(3.4) 



(3.5) 



As one can see from Fig. |2|, from the LSND region constrained by the exclusion curves of 
KARMEN and BNL-E776 and by the bound ( |3.5| ), there are four allowed regions for A^, 
Rl, R2, R3, R4 in Eq. (pj) , where we have 



A^<MU^\ 2 \U ei \ 2 <A™ x 



V4| (-', il ^ ■ (3.6) 
These regions could be explored in the near future by the MiniBooNE experiment 



Region R4 is at the limit of the final NOMAD sensitivity in the z/„ — > u e channel [56 



A mm 

|2 > Me y 



From the lower bound in Eq. Q3.6|) and the bounds (|3.1| ), we obtain 

If/,,. 12 

Therefore, we have 



+4 1 



> 



4|M 2 

A min 

4|f/ e4 | 2 



> 



A mm 
^fie 

A\TJ .12 

A mm 
A He 

4|f/e4|Lx 



(3.7) 
(3.8) 
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Since these bounds have been obtained without any assumption on the value of \U s i\ , they 
are generally valid in any 3+1 scheme. The corresponding four allowed regions for \U e 4\ 2 
and | UyA | 2 are shown, respectively, in Figs. |3| and |[ 

The lower bounds for | t/ e 4 1 2 and |^4| 2 in Eqs. ( p.7|) and (|3.8| ) imply lower bounds for the 
oscillation amplitudes B e = 4|f/ e4 | 2 (1 — |t/ e 4| 2 ) and £> M = 4|[/ M 4| 2 (1 — |£/ M 4| 2 ). Let us derive 
these bounds. 

Since |£/e4| ma x * s always smaller than 1/2, the lower bound for B e is 



mm 



A min f A 

B? n = Tfrk- i-zjTrV) ■ (3 - U) 



On the other hand, since |t^4| max can be bigger than 1/2, the lower bound for is 



BT* = min 



/tmin / /imin \ 

Me I 1 ^ e | \TT |2 (l _ ITT 12 \ 

I A 1 77 12 / ' I Mlmax V I^MImax,/ 

\ 4 l U e4| max / 



It/, 



e-1 



2 

max 



(3.12) 



Figure [| and |6| show the allowed regions in the .Be-IAmf-J and B^Am 2 ^ planes given 
by the bounds 

5f n < B a < (a = e,/i). (3.13) 

One can see that these regions lie just on the left of the Bugey+CHOOZ (Fig. |) and CDHS 
(Fig. |BD exclusion curves and could be observed in the near future [j57| . Let us emphasize 



that these regions are generally predicted in any 3+1 scheme, since they have been obtained 
independently on any assumption on the mixing (as Eq. ( |1.2| ) or Eq. ( |L3| )). 



IV. LARGE u e -> v T AND ^ -> v T SHORT-BASELINE TRANSITIONS 



We have now all the elements to calculate the bounds ( |2. 7| ) and ( |2.8| ) on the amplitudes 
of short-baseline v e — > v T and — ► z/ r oscillations that follow from the assumption ( |1.3|) , 

|f/s4| 2 <l- 

Figure ^ shows the allowed regions in the A^-IAmf^ plane given by the bounds fl2.8Q. 
One can see that the region R4 is excluded by the negative results of the CHORUS |57| and 
NOMAD |38[ experiments. The other three regions are possible and predict relatively large 



oscillation amplitudes that could be observed in the near future, especially the two regions 
R2 and R3 in which A^ T ~4x 1(T 2 - HT 1 . 

Figure || shows the allowed regions in the A^-IAm 2 ^ plane given by the bounds Q2.7|). 
One can see that these regions predict relatively large oscillation amplitudes, but unfortu- 
nately lie rather far from the CHORUS and NOMAD exclusion curves (except the region R4 
that is excluded by Fig. |7|). Therefore, even under the favorable assumption (|1.3j), it will be 
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very difficult to observe u e — > v T transitions in short baseline experiments with conventional 
neutrino beams, but large transitions could be observed with a v e beam from a neutrino 
factory ||58 



As one can see from Fig. || the bound in Eq. (|3.1|) imply that |t/ e 4| 2 is very small in all 
the three allowed regions in Fig. [7|. On the other hand, one can see from Fig. |] that |£7 M 4| 2 is 
small in the two allowed regions R2 and R3, but it is large in the region Rl. As a consequence 
of the unitarity of the mixing matrix and the assumption (|1.3[ ), we have 1 — |t/ T 4| 2 -C 1 in 
R2 and R3, whereas |t/ T 4| 2 can be as small as about 1/2 in Rl. The prediction for solar, 
atmospheric and long-baseline experiments depend on the value of [t/^l 2 . There are two 
possibilities: 

1. | VpA 1 2 "C 1 in regions R2 and R3 in Fig. In this case 1 — |f/ T 4| 2 <C 1 and the 
atmospheric neutrino anomaly is due to dominant — > v s transitions. This possibility 
is disfavored by present Super-Kamiokande and MACRO data PU| , |2"2"|| , but still it is 
not completely excluded (see |59| , |60| , |5T| ) . Obviously, in this case long-baseline — > v T 



transitions are suppressed with respect to the dominant — > v s transitions, that are 
almost entirely responsible of the disappearance of z/^'s. 

The solar neutrino problem is due to an approximately equal mixture of v e — > 
and v e — > v s transitions, which is allowed by the data. This has been shown in 
Refs. |4"6| , f4"7j in the framework of the 2+2 schemes in Fig. [TJ. However, the results 



obtained in Refs. |]46|j47j| are valid also in the 3+1 schemes, because the formalism of 
solar neutrino oscillations in 3+1 schemes is identical to that in 2+2 schemes ^-^7[. 
Indeed, from Eqs. ( |3.1| ) and ( |3.2| ) we already know that the results of the Bugey and 
CHOOZ experiment imply that |?7 e 4| 2 is very small. Furthermore, taking into account 
the hierarchy 

Am^ = Amg UN < (Am 2 ^ = Am^ < (Am^l = Amg BL , (4.1) 
the effective survival probability of v e and v e in long-baseline experiments is given by 



P^t. = 1 - MUes\ 2 (1 " |f/e3| 2 ) Sin 2 (^f^) + \U e4 \* , 



(4.2) 



where L is the propagation distance and E is the neutrino energy. Neglecting |6 r e 4| 4 , 
Eq. ([4.2|) has the same structure as the usual two-generation survival probability (see 
^]), with sin 2 2$ replaced by 4|?7 e 3| 2 (l— |f/ e 3| 2 ) and Am 2 replaced by Am^. The 



bound on sin 2 2t? obtained in the CHOOZ experiment [23] for jAm^] > 10 3 eV 2 
implies that \U e 3\ 2 < 2.6 x 10 -2 . Therefore, both \U e ^\ 2 and \U e 4\ 2 are very small, 

|t/e 3 | 2 <3x 10- 2 \U e4 \ 2 <3x 10- 2 , (4.3) 



and can be neglected in the study of solar neutrino oscillations, as done in Refs. f45|-f47 . 
In other words, the production and detection of the mass eigenstates u 3 and z/ 4 is 
negligibly small in solar neutrino experiments. Moreover, because of the hierarchy 
( f4.1| ), there are no matter-induced transitions from v\, to u 3 , Hence, i> 3 and v± 
effectively decouple in solar neutrino oscillations, leading to the formalism described 
in Refs. |5|-g7|. 
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2. (L^l — 0-33 — 0.55 in region Rl in Fig. |j. In this case the solar neutrino problem is 
due to a mixture of v e — > u^, v e — > u T and v e — > v s transitions, that is allowed by data 
as in the previous case. 

The atmospheric neutrino anomaly is due a mixture of —>■ v T and — > v s transi- 
tions. In Refs. P5| , |5H| , |5T| it has been shown that a mixture of — > v T and 



transitions in the framework of 2+2 schemes is allowed by the atmospheric neutrino 
data. This indicates that such a mixture should be allowed also in the framework of 
3+1 schemes. The groups specialized in the analysis of atmospheric neutrino data 
could check this possibility using the formalism presented in Appendix [A]. The exis- 
tence of mixed v ^ — > v T and — > v s transitions can also be checked by comparing the 
rates of z/ M disappearance and — > u T appearance in future long-baseline experiments 
with beams. The formalism that allows to describe these oscillation channels in the 
framework of 3+1 four-neutrino schemes is presented in Appendix |A]. 

These predictions are testable in future experiments, especially measuring the percent- 
age of transitions into active and sterile neutrinos in solar, atmospheric and long-baseline 
experiments (that should measure the same transitions observed in atmospheric neutrino 
experiments). Taking into account also the prediction of large — > v T and v e — > v T in 
short-baseline experiments, the schemes under consideration can be checked and possibly 
falsified in the near future. 



V. CONCLUSIONS 

In conclusion, we have considered the four- neutrino 3+1 schemes in Fig. [I], that are 
marginally allowed by present data (see Fig. 0). We have identified four small regions in the 
parameter space of neutrino mixing compatible with all data and we have derived general 
upper and lower bounds for the elements U e ± and of the mixing matrix and on the oscil- 
lation amplitudes in short-baseline v e and disappearance experiments. The corresponding 
Vfj, — > v e transitions and v e and disappearance in short-baseline experiments are relatively 
large and could be observed in future dedicated experiments. Assuming a small mixing of 
the sterile neutrino with the isolated mass eigenstate, |t/ S 4| 2 C 1, we have shown that large 

— > v T and v e — > v T transitions are predicted in short-baseline experiments. We have also 
discussed the implications of |?7 S 4| 2 1 for solar, atmospheric and long-baseline neutrino 
oscillation experiments and we have presented in Appendix [X] the formalism describing in 
general 3+1 schemes the oscillations in matter of atmospheric neutrinos and neutrinos in 
long-baseline disappearance and — > v T appearance experiments. Finally, let us remark 
that the four 3+1 schemes in Fig. [I] are equivalent for solar, atmospheric and short-baseline 
neutrino oscillation experiments, but they may be distinguishable in long-baseline — > v e 
experiments if \U e ^ is not too small ||61|1 , or through their different effects in tritium (3 decay 
and neutrinoless double-/^ decay (see @,@,H) and through neutrino oscillations in supernovae 
(see [^2|,^3j) and in the early universe (see |64H68|]). 
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APPENDIX A: FORMALISM OF ATMOSPHERIC NEUTRINO OSCILLATIONS 



In this appendix we derive in a concise way the formalism in 3+1 schemes of atmospheric 
neutrino oscillations and oscillations in long-baseline —>■ v T and disappearance exper- 
iments. We take into account matter effects, following a method similar to that used in 



Ref. [45 1 in the case of 2+2 schemes. As explained in the main text, the formalism of solar 



neutrino oscillations in 3+1 schemes is identical to that in 2+2 schemes |p5H47| . 

From Eq. ( [4.3| ) we know that and U e 4 can be neglected. Choosing the ordering u s , 
z/ e , v T for the flavor neutrino fields, the mixing matrix can be written as 

U = V U V 1A V 13 V 12 , (Al) 

where 

{Vij)ab = S a b + (cos^j - 1) {S ia 5 ib + 5 ja S jb ) + sintfy (5 ia 5 jb - 5 ja 5 ib ) (A2) 

represents a rotation by an angle #jj in the i-j plane. We have also neglected possible 
CP-violating phases. The evolution of the neutrino flavor amplitudes ip a (a = s,e, /i,r) in 
vacuum and in matter is given by the MSW equation p9[ 



i^Lq = n^, (A3) 
ax 

where \1/ = (ip s ,ip e ,ip fl ,ip T ) T and TC is the effective Hamiltonian 

H = — (U MlU ] + A) . (A4) 
2p 

Here p is the neutrino momentum, 

M 2 = diag(0, 0, Am 2 31 , Am^) (A5) 

is the mass-squared matrix, in which we neglected Am| x according to the hierarchical relation 
( (4.1| ), and A is the matrix 

A = dmg{-A NC , A cc , 0, 0) , (A6) 

with Acc = 2pVcc and A^c = ^pVnc, where Vcc and Vnc are, respectively, the matter- 
induced charged- current and neutral-current potentials (see [|2|-|6|). For atmospheric neu- 
trinos propagating in the Earth and accelerator neutrinos in long-baseline experiments 
l^ccl ~ I -A nc\ ~ l^ m 3il ^ l^w^il- in Eq. (S3) we have neglected a common phase 



for the flavor amplitudes, that is irrelevant for neutrino oscillations. 

The evolution equation ( |A"3| ) is most easily solved in the rotated basis = V^V^^/ = 
(ip'i,4 , 2,i J 3,' l Pi) T that obeys to a similar evolution equation, with 7i replaced by 

K'=^- {VisMlVg + V^AVu) (A7) 
ip 
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where we have taken into account the fact that V^-MqV 7 ^ = M.q and V34.4.V34 = A, which 
lead to a significant simplification of the evolution equation. The explicit form of 7i' is 



Ap 



( (l-c 2 ^ 13 )Am| 1 -2c| iA yl J vc s 2 tf 13 Am§ 1 s 2i } 14 A NC \ 

1 



2A C c 

i 3 Am §i ( 1 + c 2tf 13 )Am| 1 

\ -^ 14 ^JVC 2Am| 1 -2s| i4 A JV c / 



(A8) 



with c,y = cosi? and s$ = sin??. This equation shows that the amplitudes ip' 2 an d ^4 
evolve independently (remember that [Am^-J ^> |Ajvc|)j with phases given by the energy 
eigenvalues 

E' 2 = A cc /2p = Vcci ^~Am^/2p, (A9) 

whereas the evolutions of the amplitudes if)[ and ip'3 are coupled and there is a resonance in 
the 1-3 sector for 

- cos 2 $14 A NC = cos 2^13 Am^ . (A10) 

Notice that for neutrinos the left-hand side of Eq. ( |A10| ) is positive, because A NC < 0, 
whereas for antineutrinos it is negative because A^c must be replaced by A NC = —A N q- 
Hence, if Am| x > (schemes I and IV in Fig. |I|) and $13 < 7r/4 or Am| x < (schemes II 
and III) and $13 > 7r/4, there can be a resonance for neutrinos, otherwise the resonance can 
be for antineutrinos. 

One can calculate the evolution of the amplitudes ip[ and ip' 3 solving numerically the 
two coupled equations generated by the 1-3 sector of the effective Hamiltonian 7i' . Another 
common method for the solution of the evolution equation is to divide the Earth interior 
into shell with constant density, calculate the evolution of the amplitudes in each shell 
and match the amplitudes in the flavor basis at the shell boundaries. In this case, in each 
shell the amplitudes ip[ and i/j 3 can be written as 

^ = cos tffgVf + sin , tp' 3 = - sin tf^Vf + cos <^ 3 M , (All) 

where ipi 1 an d ^3 are the amplitudes of the energy eigenstates that evolve with phases 
given by the energy eigenvalues 



P M _ 1 
£j\ 3 — — 

M Ap 



Ami - cos 2 $14,4 jvc =f \l (cos 2$ i3 Am^ + cos 2 $u A NC f + (sin 2$i 3 Am? ^ 2 



(A12) 

The effective mixing angle in matter is given by 

„w sin 2i?nAmo, 

tan 2^ 3 = o n . „ 13 — - , A13 

13 cos 2$ 13 Aml 1 + cos 2 14 A NC v ' 

from which one can see that "&^ 3 = tt/A when the resonance condition ([AlOj) is satis- 
fied, and ^ — > 7r/2 in a very dense medium, where | cos 2 ^uA NC \ 3> | cos2t? 1 3Am3 1 | if 
cos 2$i3 Am 2 ! > 0, whereas tfjf — > in a very dense medium if cos2$i3Am 2 1 < 0. 
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The connection between the amplitudes ty' and the flavor amplitudes needed for the 
calculation of the probability of flavor transitions, is given by 



= cos #14^4 + sin $i4^/>4 , 



ipu = — sin $ 



. j i — n fluip'i + cos i?34'03 + cos i?i4 sin $34t/> 4 , 

= — Sin $14 COS 1?34'0i — Sin $34^3 + COS $14 COS $34^4 • 



(A14) 



The second line of Eq. (|A14|) imply that electron neutrinos do not oscillate in atmospheric 
neutrino experiments (remember that ip'2 evolves independently). This is due to the approx- 
imation U e s = 0, motivated by Eq. ( |4.3|) , which also imply that the charged-current matter 
potential Vcc-, felt only by u e , is irrelevant for atmospheric neutrino oscillations (only the 
neutral-current potential Vnc, felt by u n an d v T , enter in Eqs. (|A10|) - (|Al3|) ). On the other 



hand, simultaneous — > v T and — > v s transitions are allowed, with pure — > iv transi- 
tions in the limit cos$i 4 = 0, that corresponds to [7 s4 = 1, and pure — > za, transitions in 
the limit cos $ i4 cos $ 34 = 1, that corresponds to {7 T 4 = 1. 

Let us consider finally long baseline experiments in which the neutrino beam travels in 
the crust of the Earth, where the matter density is practically constant. The probability of 
v a transitions is given by 



LBL 



+ "/3 



U al U f3l 



TjMjjM 
U a3 U 03 eX P 



-ATJ M TJ M TJ M TJ M sin 2 
q:U al U 0l U a3 U 03 Sm 



2p 



AffL 
Ap 



M 1 2 1 t jM 1 2 



TJM VI TT 

u a4 \ I U , 



04 \ 



rjl ttM \2\ttM\2 
Z \ U a4\ \ U 04\ 



(A15) 



where L is the propagation distance, 



= 2p (£ 3 M - El 1 ) = V /(cos2$i3Am|i+cos 2 $i 4 ^ c ) 2 + (sin2$i 3 Am|i) 2 , 



(A16) 



and U is the effective mixing matrix in matter, 



U 



M 



( 



\ 



^14 S -&™ 



\ 



»M-^,.C M C,, C# 



/ 



(A17) 



where is equal to Vi 3 with $ 13 replaced by $ff . The expression (|A15[ ) can be used to 
analyze the data of long-baseline z/^ — >■ z/ T and disappearance experiments (the analysis 
of long-baseline — *> v e data requires the relaxation of the approximation U e % = 0, leading 



to a significant complication of the formalism [51 1) 
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FIG. 1. Qualitative illustration of the possible four- neutrino schemes. 
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FIG. 2. Very Thick Solid Line: Allowed regions. Thick Solid Line: Disappearance bound 
Dotted Line: LSND 2000 allowed regions at 90% CL gl. Solid Line: LSND 2000 allowed 



regions at 99% CL pi]] . Broken Dash-Dotted Line: Bugey exclusion curve at 90% CL [32]. 
Vertical Dash-Dotted Line: CHOOZ exclusion curve at 90% CL Long-Dashed Line: 

KARMEN 2000 exclusion curve at 90% CL pi Short-Dashed Line: BNL-E776 exclusion curve 
at 90% CL iF 
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FIG. 3. Dotted Line: \U e 



obtained from the 90% CL exclusion curve of the Bugey 



short-baseline reactor v e disappearance experiment |32| (Eq. ( |3,2|) ), Dashed Line: |C/ R 4| 2 nav ob- 



'e4|max 

; e4|max 

tained from the 90% CL exclusion curve of the CHOOZ long-baseline reactor v e disappearance 
experiment pj| (Eq. (|3^)). Solid Line: Allowed regions (Eq. (|3.9|) ). 
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obtained from the exclusion curve of the CDHS accelerator 



disappearance experiment Q (Eq. ([T^)). Dashed Line: |£^4| max obtained from the up-down 
asymmetry of multi-GeV muon-like events measured in the Super-Kamiokande experiment |2 ( 
(Eq. (U)). Solid Line: Allowed regions (Eq. (glcD ). 
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FIG. 5. Solid Line: Allowed regions. Dotted Line: Bugey exclusion curve at 90% CL [32] 
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FIG. 6. Solid Line: Allowed regions. Dotted Line: CDHS exclusion curve at 90% CL [33]. 
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FIG. 7. Solid Line: Allowed regions. Long Dashed Line: CHORUS exclusion curve at 



90% CL Short Dashed Line: NOMAD exclusion curve at 90% CL @. Dotted Line: 

CDHS exclusion curve at 90% CL 
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FIG. 8. Solid Line: Allowed regions. Long Dashed Line: CHORUS exclusion curve at 



90% CL [f37j]. Short Dashed Line: NOMAD exclusion curve at 90% CL f3§. Dotted Line: 



Bugey exclusion curve at 90% CL [32|. Dashed Line: CHOOZ exclusion curve at 90% CL |23 ]. 
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